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abstract : The possible sources of light nuclei populations observed recently below the geomag- 
netic cutoff by the AMS experiment are discussed in terms of nuclear processes: fragmentation of 
the incoming flux of cosmic 4 He on atmospheric nuclei, and nuclear coalescence from proton and 
t— I ' A He induced reactions. Results of simulations for 2 ' S H and 3 ' 4 ife, are presented and compared 

to the data. 

CN ' The study of particle populations in the earth environment has a long and rich history covering the 
J>y last few decades (see [1] for a general overiew of the subject, and the references in [2, 3] for details). 
After a period of waning activity, the topics is likely to regain interest with the occurence of a new 

^ ' generation of balloon and satellite experiments, wich open prospects for data samples of unmatched 
quantity and quality. This is well illustrated with the new data obtained recently in orbit close to earth 
by the precursor flight of the AMS experiment, which achieved on its orbit, a set of high accuracy 
measurements of charged particle flux over a wide range of latitude. New features of the proton [2, 4] 

■ and lepton [3] flux have been uncovered, and surprisingly large values of flux measured below the 
geomagnetic cutoff (GC). These features have been successfully interpreted in previous works by the 

| authors [5, 6], referred to as I in the following. In addition, small but significant populations of 3 £feand 

■ deuterium (D or 2 H) particles were also measured below GC, with kinetic energies extending beyond 
1 GeV per nucleon [7, 8]. This paper completes a set of three reports devoted to the interpretation of 
the new AMS data. Its purpose is to investigate the possible origins of the measured 3 He and D flux 

■ using the same approach as reported in I. 

Since particles below GC cannot be primary cosmic rays (CR), they have to be produced by nuclear 
reactions between incoming CRs (mainly p and A He), and atmospheric nuclei (mainly 14 N and 16 (9). 
The pattern of Z=2 particle spectra observed above and below GC is highly peculiar since only 4 He 
axe observed above GC whereas only 3 ffeare found below GC, with however a small admixture of 
the other isotope being compatible with the data in both cases [7]. This pattern, together with the 
. £h ! relative population of light nuclei, provide clues to the dynamical origin of the subGC particles (or 
^ ■ Albedo particles in the geophysical terminology). It can be understood qualitatively and evaluated 
quantitatively in terms of the nuclear reaction mechanisms involved in the production process. 

Collisions between nuclear systems such as those of interest here, for incident energies beyond the 
nuclear fermi energy (~ 35 MeV/nucleon) region, have very characteristic distributions. The rapidity 
distributions of reaction products at forward angles display two peaks centered around projectile 
and target rapidities, corresponding to projectile and target fragmentation in peripheral collisions 
respectively, and an intermediate plateau between these limits stemming from more central collisions 
(see for example [9]). Experimentally, the (projectile) fragmentation regime sets in as low as 20 MeV 
kinetic energy per nucleon [10]. The features of the measured differential cross sections at small angles 
can be accounted for in the fragmentation model [11], with the width of the fragmentation peak being 
governed by nuclear motion and nucleon arrangement probability. 

With the increasing production angles, the target-like and projectile-like peaks move to the mid- 
rapidity region, leaving a single broad mid-rapidity peak surviving at the largest angles [9]. This latter 
kinematical domain is associated to smaller impact parameters and larger energy dissipation due to 
larger density overlaps in the collisions. In this region, the simple fragmentation picture [11] fails and 
the description of the collision in terms of the various models based on thermodynamics and spectator- 
participant pictures [12, 13, 14] are more appropriate The energy distributions of very light fragments 
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like 1,2,3 H, 3,4 He, produced at large scattering angles in such collisions, can be described by means of a 
variety of models which are all variations approximating more or less successfully the complex multiple 
nucleon-nucleon scattering processes and nuclear collective effects governing these collisions [12, 15, 16] 
(see also [13, 14, 17]). In this context, the coalescence model [18, 19] stands as the most successful 
phenomenological approach of the production cross section for light nuclei, with a remarkable ability 
to reproducing data over a very wide range of incident energies [14, 17, 20, 21]. In this approach, 
nucleons coalesce into clusters whenever they fall within the coalescence radius (momentum) in the 
final state of the collision. Note that the models used here are considered from a purely practical 
point of view, and not be discussed in their foundations. 

For the present purpose, the phenomenological perspective can be summarized the following way: 
Small production angles are dominated by velocity-conserving projectile-like fragmentation products, 
while at large angles the coalescence production of light fragments of much lower energy per nucleon, 
dominates the cross section. In the fragmentation picture, the cross-section is expected to decrease 
with the decreasing mass number of the (projectile -like) fragment (larger probability for smaller 
number of nucleon transferred [10, 22]), whereas conversely, in the coalescence model it is decreasing 
with the increasing fragment mass number (larger probability for smaller cluster mass) . In the case of 
4 He projectile, projectile-like particles are also light fragments, then likely to be produced either by 
fragmentation or by coalescence, wheareas obviously, protons can only induce production of coalescence 
fragments. 

Projectile fragmentation in heavy ion collisions has been extensively studied experimentally in the 
past [23]. A few 4 He fragmentation studies have been reported in the literature over the momentum 
range of interest here [24, 25, 26]. The production of non velocity-conserving light fragments in p and 
4 He induced collisions is also fairly well documented experimentally [13, 17, 14]. 

The s He spectrum observed in the AMS experiment [7] between the low energy cutoff of the spec- 
trometer and up to about 1 GeV/nucleon (see figures 2 and 3 below), could originate from the two 
types of reactions described above: 

1) It could be produced by CR 4 He particle fragmentation or fragmentation-like process on atmospheric 
nuclei (this assumption was also quoted recently in [27]). The quasi absence of subGC 4 He is easy to 
understand in this framework since secondary 4 He from fragmentation could only originate from the 
very small CR flux of heavier nuclei, i.e., mostly 12 C and 16 (the possible 4 He yield from nuclear 
evaporation being expected at energies below the AMS sensitivity range). However, the basic velocity- 
conserving property discussed above, is clearly not met here since incident CR 4 He have momenta 
larger than about 6 GeV per nucleon in the equatorial region, while detected 3 He and D particles 
have less than 1 GeV per nucleon (see details below). 

2) The observed s He population could also be produced by nuclear coalescence [18] (see [14] for a review 
of models) from p and 4 He induced collisions on atmospheric nuclei [17, 28, 29]. It would explain as 
well the absence of subGC A He as discussed above, with a larger population expected below GC for 
3 ^fethan for 4 He particles, by typically one order of magnitude [17], which is compatible with the 
experimental observation. 

The experimental situation then favors the coalescence assumption. The measured flux however, 
depend on the particular dynamics of each production reaction and of the acceptance of the magne- 
tosphere to the particle considered, and only a detailed investigation can provide a definite answer. 

The inclusive spectrum of light nuclei flux at the altitude of AMS (390-400km) have been calcu- 
lated by means of an evolved version of the same simulation program as described in I. CR protons 
and Helium 4 are generated with their natural abundance and momentum distributions. They are 
propagated inside the earth magnetic field using 4th order adaptative Runge-Kutta integration of the 
equation of motion. They are allowed to interact with atmospheric nuclei and to produce secondary 
nucleons p,n, 2,3 H, and 3,4 He, particles with cross sections and multiplicities as discussed below. Each 
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secondary particle is then propagated and allowed to interact as in the previous step. Only destructive 
interaction is taken into account for light nuclei, except 4 He . A reaction cascade can thus develop 
through the atmosphere. The reaction products are counted when they cross the virtual sphere at the 
altitude of the AMS spectrometer, upward and downward. All charged particles undergo energy loss 
by ionisation. Each event is propagated until the particle disappears by nuclear collision, stopping 
in the atmosphere by energy loss, or escaping to outer space. See I for details. The small produc- 
tion cross section, and then multiplicity, combined with the small magnetosphere acceptance, for the 
particles of interest here, require a tremendous number of events to be generated, and then a huge 
computer time, for significant statistics to be reached. This difficulty has been turned around by 
enhancing numerically the production multiplicity, with the produced events being weighted by the 
inverse enhancement factor. It has been checked carefully that no distortion of the physics involved 
could be induced by using this method, in the studied case. 

The CR proton and helium flux used in the calculations were those measured by AMS [2, 4, 7]. 
The values of the total p and 4 He reaction cross sections used were based on the parametrization of 
[30], the latter being checked on the measurements performed on carbon target [31]. The production 
cross sections for light nuclei have been implemented and run simultaneously in the event generator 
on the basis of the two models described above. 

A) - For the 4 He fragmentation cross section, the model from [11] was used. In this model, the 
fragment production cross section is proportional to e in the projectile reference frame, with 

9 9 A f ( A -Af) 

a = Oq a p _ 1 — , and P, A p , Af being the fragment momentum in this frame, projectile mass number, 
and fragment mass, respectively. The value of Fermi momentum related parameter o"o=100 MeV/c 
was chosen in agreement with the results of [11, 22]. The total A He fragmentation cross-sections into 
3 He , 3 H, and D used in the program, were borrowed from measurements on Carbon target [32], and 
corrected for their A 1 / 3 dependence. 

B) - In the coalescence model, the invariant differential production cross sections for composite frag- 
ments with mass A, are related to the nucleon production cross section by a simple power law: 

(F / 3 N A d 3 N p A 

dp A dp p 

with p~A = A ■ p~p. This relation provides straightforwardly the momentum spectrum of mass A 
fragments as a function of the nucleon spectrum at the same production angle. The inputs of the 
event cross section calculation then consist only of the value of the Ba parameter and of the proton 
production differential cross-section. The coalescence parameters for p induced collisions have been 
found to have approximately constant values through the energy range from 0.2 GeV/n up to 70 GeV/n 
[33]. The following values, averaged from [14, 17, 21, 29, 34], were used in the program: £>2 = 2.5- 10~ 2 
for deuteron, B 3 = 2.5 • 10" 4 for 3 He , B' 3 = 4 ■ 10~ 4 for 3 H, and B A = 4 • 10~ 6 for 4 He , in units 
(GeV 2 /c 3 )( A_1 ). The accuracy on these values is estimated to be within ±30% for £>2, and ±50% for 
B% and B 3 . See [17] for B A . The proton induced proton production was generated as described in [35] 
using the paramertization from [36], while for incident 4 He particles, the same spectral shape as for 
protons was used, with the appropriate scaling to take into account the experimental total reaction 
cross section [31] and proton production cross section (multiplicity) for A He collisions on nuclei [26, 37]. 

The simulation has been run for 2 10 8 incident primary proton and helium particles generated 
at the injection sphere. This number corresponds to a sampling time of 4 10~ 12 s of the cosmic ray 
flux. The results are shown on figures 1 to 4. Note that no adjustable parameter was involved in the 
calculations. 

The general features of the simulated sample are quite similar to those already reported in I for 
protons. In particular, concerning the trapping of particles in the earth magnetic field, the same trend 
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is observed for the population to be confined with a long lifetime (> 10 s) in the region of equatorial 
latitudes. The same class of low energy very long lifetime but low crossing multiplicity, quasi-polar 
population as observed in I, is found here, likely corresponding to the population of more outer belts 
than studied here. 

Figure 1 shows the some basic distributions of the produced 3 ffeand D particles, crossing the 
detection altitude (no angular acceptance involved at this level). The fractional energy of the produced 
fragment with respect to the energy of the incident article (top left) shows a clear distinction between 
the coalescence yield at low energy and the approximately velocity-conserving fragmentation products 
at high energy. Also seen is the strong dominance of the D coalescence yield, larger than the 3 He yield 
by two orders of magnitude, and the much less different fragmentation yields, due to the particular 
structure of 4 He which fragments into D with a naturally higher multiplicity than in the general 
case. The rank distribution in the atmospheric cascade for 3 He (top right, solid histogram) displays 
a prominent peak for rank one. This originates from the fragmentation 3 He component which occurs 
mostly at the first interaction of the primary CR 4 He . Coalescence particles are seen to be produced 
up to more than 8 generations of collision in the cascade. The mean production altitude (bottom 
left) predicted for deuterons (40 km) in agreement with our previous calculations (see [38] for the 
experimental context) is significantly lower than for (coalescence plus fragmentation) 3 He(51 km). 
This is also due to the fragmentation process which occurs at the first collision and thus at higher 
altitude on the average (dotted histogram). The 3 He coalescence yield (not shown) has a similar 
distribution as for Ds. The same strong East- West effect (bottom right) responsible for the lepton 
asymmetry reported in [6], is also seen here as it could be expected. 

Figure 2 shows the comparison between the flux spectra measured by AMS (full circle) [7] and 
the simulation results (histograms), taking into account the spectrometer acceptance (30 deg), for 
3 bins in latitude. The solid histograms correspond to the 4 He CR flux above GC. It is seen that 
they reproduce fairly well the experimental distributions, with however a tendancy to underestimate 
the experimental CR flux close to GC. This defect was not observed in our previous works on the 
protons [5] and lepton [6] flux. The dotted histograms correspond to 3 He particles produced by 4 He 
fragmentation. The expected yield is seen to be significant only above GC, and the differential flux 
appears to be more than two order of magnitude smaller than the primary 4 He flux, except for the few 
bins very close to the cutoff. This value is small compared to the known 3 He CR flux [39, 40, 41] which 
ratio to the CR He flux is about 10% for this momentum range, i.e. about ten times larger than the 
value calculated here. This result is compatible with the AMS measurements (see figure 4 in [7]). The 
dashed and gray-shaded histograms correspond to CR p plus A He induced coalescence 3 He and 4 He 
flux respectively. As expected from the introductory discussion, the 3 He flux is more than one order of 
magnitude larger than for 4 He . These flux are lying below GC, and the simulated 3 He spectra account 
pretty well for the measured spectra in magnitude and shape, with some underestimate of the data 
however in the higher subGC momentum range for the equatorial and intermediate latitudes. 

Upper figure 3 shows the (downwards) 3 He particle spectrum measured by AMS, in kinetic energy 
per nucleon, integrated over the corrected geomagnetic (CGM, see I for the precise definition) latitude 
latitude range \6m\ < 0.6 rad, compared to the simulation results (coalescence yield, solid histogram). 
The calculated values are seen to be in agreement with the data within a factor of about 2 over the 
whole energy domain and over a dynamical range covering two orders of magnitude, which can be 
considered as a very good overall agreement. The small 4 He coalescence yield (dashed histogram) 
predicted is expected to be more than one order of magnitude smaller, a value compatible with the 
AMS conclusions in which an experimental upper limit of 10% for this ratio was set. Lower figure 3 
compares the experimental (full circles) and simulated (histograms) energy-integrated 3 He flux below 
GC as a function of the geomagnetic latitude. The full histogram corresponds to the total expected 
yield. It accounts quite well for the measured flux except for the highest latitude bin. The dashed 
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histogram gives the 4 He induced contribution, which is seen to amount to 40-50% of the total yield. 
This is due to the combination of a larger total reaction cross section and larger proton multiplicity 
(with then a larger coalescence probability) with 4 He than with proton incident particles, which 
enhances the 4 He induced coalescence yield. 

Figure 4 shows the expected spectral yields for 2 H and S H particles, for a set of bins in latitude. 
Note that although the two contributions were included, these yields are almost exclusively due to 
coalescence. The results for the 2 H spectra below GC are in very good agreement in magnitude and 
in shape with the preliminary AMS results [8]. This ability of the calculations to successfully and 
consistently reproducing both the 3 He and D flux at the same altitude of measurement with the same 
production model constitutes a very strong indication that these particles do originate from the same 
coalescence mechanism. 

Some 3 He flux measurements have been reported recently at lower energies and larger distance 
from earth in the region of the inner belt [42]. This flux could not be definitely interpreted in the 
quoted work. It could certainly be investigated in the present approach. 

In conclusion, it has been shown that the AMS measurements of the s He and D particle flux can 
be reproduced consistently and simultaneously, together with the proton secondary flux as reported 
in I, by a simulation incorporating the interactions between Cosmic Ray flux, earth magnetosphere 
and atmosphere, and assuming the 3 He and D particles to be produced by coalescence of nucleons in 
the CR proton and 4 He induced nuclear collisions with atmospheric nuclei. The 4 He fragmentation 
products appear not to contribute at a detectable level to the flux measured below the geomagnetic 
cutoff. 

These results also clearly point to the interest for the future satellite experiments to have a particle 
identification capability covering this mass region and extending over a large kinetic energy range down 
to about 0.2 GeV per nucleon, in order to allow the experiments to collect data samples larger by more 
than two orders of magnitudes than those analyzed here, and then to make possible the achievement 
of a much more detailed investigation of the issues addressed in the present work. 
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Figure 1: Features of the simulated samples for 3 He (dashed histograms) and D particles (solid histograms). 
From top to bottom and from left to right: Altitude of production (km), rank in the atmospheric cascade, 
fractional energy of particles, and Est-West angle distributions (rad). 
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Figure 2: Comparison of the Helium spectra measured by AMS ([7] full circles) at various CGM latitudes 
9m, with the simulation results (histograms). Full line : He above GC; Dotted: 3 He from He fragmentation; 
Dashed: coalescence He ; shaded: coalescence ^He . See text for details. 



9 



> 

CD 
O 

1/2 

O 
CD 



10 



A 10 



!0J < 0.6 



1 



10 



0.02 



0.2 0.4 0.6 0.8 1 1.2 
Kinetic Energy per Nucleon (GeV/A) 



Sh 

CJ 
CD 

1/3 



X 



0.01 



i 1 1 1 1 1 1 

0.1 0.2 0.3 0.4 0.5 0.6 
Corrected Geomagnetic Latitude 

Figure 3: Top: Experimental energy spectrum of the subGC He flux integrated over 9m < 0.6 rad, measured by 
AMS [7], compared to simulation results (solid histogram) . The dashed histogram corresponds to the coalescence 
4 He flux. Bottom: Energy-integrated 3 He experimental distribution (E < 1.2 GeV/A) as a function of the CGM 
latitude (full circles) compared to simulation results for proton and A He induced flux. The dashed histogram 
shows the He induced contribution. 
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Figure 4: Prediction for the subGC flux of deuterium 2 H (thick histograms) and tritium 3 H particles (thin 
histograms) for various bins of latitude between equator and polar region (CGM bin values 9m given in radian). 
Solid line : full (p plus 4 He ) yield; Dashed line : 4 He yield. 
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